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Abstract 



We study a low-energy effective field theory (EFT) describing the NN system 
in which all exchanged particles are integrated out. We show that fitting 
the residue of the 3 Si amplitude at the deuteron pole, rather than the 3 Si 
effective range, dramatically improves the convergence of deuteron observables 
in this theory. Reproducing the residue ensures that the tail of the deuteron 
wave function, which is directly related to NN scattering data via analytic 
continuation, is correctly reproduced in the EFT at next-to-leading order. The 
role of multi-nucleon-electroweak operators which produce deviations from 
effective-range theory can then be explicitly separated from the physics of the 
wave function tail. Such an operator contributes to the deuteron quadrupole 
moment, (j,q, at low order, indicating a sensitivity to short-distance physics. 
This is consistent with the failure of impulse-approximation calculations in 
NN potential models to reproduce hq. The convergence of NN phase shifts 
in the EFT is unimpaired by the use of this new expansion. 



ust 1999 



When the deuteron is probed at low momentum-transfers it is appropriate to treat the 
nuclear bound state in a long-wavelength approximation. In such an approximation no 
details of the short-distance dynamics are observed. Consequently most observables are 
dominated by the physics of the deuteron's wave function tail. Ignoring, for the time being, 
the deuteron's D-state, two numbers characterize this tail. One is the deuteron binding 
energy, which determines the rate of exponential fall-off of the wave function. The other is 
the asymptotic S-state normalization, As, which multiplies this exponential. In the limit as 
r —>■ oo the wave function becomes 



^ 1 r V 2tt r K } 

Below we work in terms of the dimensionless quantity Zd, rather than As = y/2r/Zd. 
Neither Zd nor As can be measured directly in experiment, but their values can be deter- 
mined by analytic continuation of the NN scattering amplitude into the bound-state region. 
Parameterizing the 3 Si phase shifts as JIJ 

|k|cot5 = -7 + ^(|k| 2 + 7 2 ) + W2 (|k| 2 + 7 2 ) 2 + - , (2) 

with 7 _1 = 4.319 fm, pd = 1.765(4) fm, and w<i = 0.389 fm 3 0[§, reproduces the phase- 
shift data in the low-momentum region. Near the deuteron pole, the scattering amplitude 
becomes 
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Here the quantity R(k) is regular in the limit k — > 27 and, importantly, is the only place 
that higher-order terms in the expansion of |k| cot^o enter this amplitude. It follows that 
the residue Z d is given by 
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The measured values of 7 and pd quoted above give Z d = 1.690(3), which is also the value 
of Z d obtained using A s given in the 1993 Nijmegen phase-shift analysis [0|§, as required. 

In the 50 years since effective range theory was first formulated many calculations have 
taken the "effective-range-theory wave function" given in eq. (|1|) as the wave function for all 
r and calculated low-momentum processes involving the deuteron. As we will recapitulate 
below, the results obtained by this procedure are often in good agreement with both exper- 
imental data and modern NN potential models. This suggests that such quantities truly 
are insensitive to details of the short-distance deuteron dynamics. 

The central tenet of effective field theory (EFT) is that low-momentum probes of any 
system should not reveal details of that system's short-distance dynamics. Much effort has 
recently been invested in calculating observables in the two-nucleon system using effective 
field theory techniques based on this tenet (§]- An extreme, but efficacious, EFT has 
been constructed for the NN system by considering only distance scales r l/m n , i.e. 
the region where the S-state deuteron wave function is given by eq. (H). In this EFT the 



2 



low-momentum scales in the problem are 7 and the typical nucleon momentum k. It is a 
non-relativistic theory of nucleons interacting via contact interactions, which we denote by 
EFT(^"). These contact interactions represent the results of integrating the pions, and all 
other exchanged degrees of freedom, out of the theory. The leading contact interaction is 
responsible for binding the deuteron, and it is enhanced over expectations based on naive 
dimensional analysis, as are all other S-wave contact interactions. Relativistic corrections 
can be included perturbatively, and are very small. We shall not discuss them here, the 
interested reader should consult Ref. |I6| for details. The Lagrangian describing the strong 



interactions in EFT (7^ is 

£ = -ifcf^ (N T P t N) j (N T P t N) - *ct Sl) \ [(N T P t Ny (N T O^N) + h.c] + ... , (5) 



where Pi is the spin-isospin projector for the 3 Si channel [|13|,[nj and the ellipses denote 
terms involving more spatial derivatives. The two derivative operator, O^' 3 , is given by 
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The coefficients Co, C%, ... are to be fit to data. Until now this has been done by repro- 
ducing the effective-range expansion (^) up to some given order in the expansion parameter 
Q ~ k, 7. This is straightforward, since the exact result for |k| cot 5 in EFT(^"), if relativity 
and S-D mixing are ignored, is 

47T 1 

- |k|cot(5o = +/i , (7) 



when dimensional regularization with power-law divergence subtraction is employed JIB" 
Expanding each C<i n in powers of Q itself, i.e. writing: Cm = Gm-n-\ + Cm-n + • • •, and 
then expanding eq. (0) in powers of Q one can match to eq. (^) order- by-order in Q. If pd 
and W2 are both taken to scale as Q°, then at next-to-next-to-leading order (N 2 LO) in Q 
this yields the coefficients 

f r ( 3 Si) _ _ 4tt 1 i (B Sl) _ 2vr j 2 p d . ( 3 Sl) 7T 7 4 P1 

" M N - 7) ' " M N 0* - 7 ) 2 ' ^ " M N Qi- 7 ) 3 ' 

^ 2 ~ M N (p - 7 ) 2 ' ^ 2 - 1 ~ M N (p- 7 ) 3 ' A '~ 3 " M N (p- 7 ) 3 ' (8) 

at the renormalization scale p. The dominant interaction ^Cq f 1 ^ must be iterated to all 
orders to generate the leading order (LO) amplitude. All other interactions, including the 
corrections to Co, can be treated in perturbation theory. At next-to-leading order (NLO) 
and above the coefficients are chosen to reproduce pd exactly. Hence we refer to this choice 
as the "p-parameterization. With this choice of C's, observables are power-series expansions 
in the quantity •ypd order-by-order in the Q expansion. In particular, the expansion of Z d is 

Zf ] = l +1Pd + ( 1Pd f + ( 1Pd f + ... . (9) 

Physically, the amplitude for any low-energy elastic reaction on the deuteron will always 
include an overall factor of Zd- Thus the appearance of this series in a number of previous 
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calculations in EFT (ft) should come as no surprise. Indeed, Park et al. have had con- 
siderable success in reproducing experimental data by pursuing a strategy of fitting Z d at 
NLO 0,0. In fact, such an approach is implicit in other work on EFT in nuclear physics, 
e.g. Refs. [^,10|)29|. Recently, Phillips and Cohen |4l|] and Rho [4^] have stressed that it is 



better to fit Zd in the EFT than pd, since this ensures that the calculation correctly repro- 
duces the long-distance piece of the deuteron wave function. The difference between fitting 
Pd and Z d is higher-order in any given EFT calculation, but, as we shall show here, in EFT^) 
fitting Zd markedly improves the convergence of calculations for a variety of low-energy elas- 
tic processes on deuteron targets. We will also discuss inelastic processes such as np — > dj, 
and show that demanding that Zd is correctly reproduced also improves the convergence of 
these calculations. Finally, we will discuss the NN phase shifts in the Q-expansion when 
this alternative fitting procedure is employed. We will show that, despite abandoning the 
exact fit of pd in the Q-expansion of k cot S, the results found for the phase shift with this 
modified fitting procedure are no worse than those found using the "p-parameterization" . 

In this work we require that the position and residue of the deuteron pole are reproduced 
exactly at NLO in the Q expansion. The coefficients of the Lagrangian ([5]) determined by 
this constraint are (again neglecting relativistic effects): 

i r esi) _ 4?r 1 # r ( 3 Si) _ 2?r -i(Zd-l) t r (3 5l) _ 7T (Z d - l) 2 (2p- 7 ) 7 

" M N (p - T ) ' ^°>° " M N (p- 7 ) 2 ' ^ ~ M N (p- 7 ) 3 

^( 3 Si) - 2n Zd ~ l ir^ _ 2n ( Zd ~ X ) 2/i ir^ 71 ( Zd ~ l ^ (m\ 

Z 2 - 2 ~M Nl (p-^ ' ^ " M^ 7 (^-7) 3 ' " M^7 2 (^-7) 3 ' ( } 

We call this the "z-parameterization"[|. Here the quantity Zd — 1 is taken to be of order 
Q. This choice for the C's is formally equivalent order-by-order to that of eq. @, the 
difference between the two always being of higher order in the Q-expansion. The advantage 
of this choice can be seen in the Q-expansion of Zd- Instead of the convergence displayed in 
eq. (|9[) we now have, by explicit construction: 

zf ] = l + (Z d -l) + + + + ... . (11) 

Both parameterizations yield Z d L °^ = 1. However, at NLO the two series in eqs. (||) and 
(|TT|) have very different behavior. In the p-parameterization the first 7p^ correction brings 
that expansion perturbatively close to the complete sum of Zd, with corrections of order 
(lpd) 2 and higher. In contrast, in eq. (|TT]) the first (Zd — 1) correction brings the series into 
exact agreement with the complete sum, and all higher-order terms vanish. This behavior 
translates directly into the convergence of the Q-expansion for observables involving the 
deuteron. 

The easiest place to see this is in the deuteron charge form factor Fc- In the p- 
parameterization, neglecting relativistic effects gives LO, NLO, and N 2 LO contributions 
to the charge form factor which are [ftB], 



Pronounced "zed" -parameterization . 



4 



F c°\H) = ^n~ 1 (JST) , F^(|q|) = - 1Pd 
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where |q| is the momentum transfer, and = 0.790±0.012 fm is the isoscalar charge radius 
of the nucleon. In contrast, the ^-parameterization yields LO, NLO and N 2 LO contributions 
to the form factor: 





^ (2) (|q|) = -^|q| 2 ^tan-M^] . (13) 



Comparing the expressions for Fc in the p and z parameterizations it is clear that 
demanding that be reproduced exactly at NLO means that the only correction at N 2 LO 
is due to the finite size of the nucleon. There are no corrections resulting from "wave 
function" effects at any order beyond NLO. 

This can be clearly understood if we assume the wave function given in eq. ([!]), is correct 
for all r, calculate Fc via: 

F c (|q|)= [ d 3 re^ 2 \ij(r)\ 2 , (14) 



and also add an additional (constant) piece to Fc which ensures that Fc(0) = 1. When this is 
done we reproduce eq. (|IB|), but with r N = 0. In other words, at NLO the ^-parameterization 
exactly reproduces the results obtained in effective range theory calculations of a deuteron 
containing point-like nucleons. However, EFT also allows for systematic incorporation of 
the effects of finite-nucleon size, via terms such as those proportional to r 2 N above, in a way 
that maintains the locality of the field theory. 

Furthermore, one order beyond that calculated here, i.e. at N 3 LO, a four-nucleon-one- 
photon interaction contributes to Fc and shifts the deuteron charge radius from its impulse 
approximation value. This contribution to Fc is beyond the scope of effective range theory. 
In contrast, the charge radius extracted from eq. (0) at N 2 LO is 

(rT FT = r 2 N + #| = 4.566 fm 2 . (15) 

87^ 

We estimate the error from the next-order corrections (including the counterterm) to be 
~ 3%. In Table | we compare this, and other static properties of the deuteron, to those 
obtained in the potential models Nijm93, Reid93 [fPJ, and OBEPQ |H5f (energy-independent 



version), in impulse approximation, and also to the experimental value. We see that in 
the case of (r 2 .) the EFT calculation agrees with data and potential model calculations 
within 1-2 %. Note that in comparing results we have assumed a common contribution of 
r 2 N = (0.79 fm) 2 to (r 2 d ). 

In Table H we have also displayed results for the deuteron's quadrupole moment, /xq, and 
for the quadrupole "radius", (tq). We now spend some time discussing the results obtained 



5 



Calculation 


(rg) (fm 2 ) 


m (fm 2 ) 


{rl) (fm 4 ) 


EFT (LO) 


2.332 


0.335 


0.70 


EFT (NLO) 


4.566 


0.286 (fit) 


1.40 


OBEPQ 


4.499 


0.274 


1.27 


Nijm93 


4.489 


0.271 


1.22 


Reid93 


4.501 


0.270 


1.22 


Experiment 


4.493(23) 


0.2859(3) 





TABLE I. Results for the square of the deuteron charge radius, the deuteron quadrupole 
moment, [1q, and the square of the quadrupole "radius", {tq), at LO and NLO, compared with im- 
pulse-approximation potential-model calculations, and the experimental results. The experimental 
value of (r?) is taken from Ref. [46|. See also the more recent Refs. p7|,E3]. 



in EFT (ft) for the deuteron's quadrupole form factor. This quantity clearly demonstrates 
the limitations of effective range theory. 

Naively one might think that it is impossible for the deuteron to acquire a quadrupole 
moment in this EFT, but such is not the case. The deuteron becomes non-spherical due 
to the appearance in the theory of contact operators with two derivatives which mix the S 
and D-states. In Ref. it was shown that the coefficient of this operator could be fit to 
the asymptotic D-to-S ratio, r] s d. The value of r] s d is known to high precision from analytic 
continuation of the S-D mixing parameter, ~E\\ r\ s d = 0.02544 ± 0.0002 [|J. Fitting T] S d means 
that the tail of the D-wave piece of the deuteron wave function will be correct. The leading 
order result for Fq is then just the same as that obtained assuming asymptotic forms for 
the deuteron radial wave functions, working to leading order in rj s d, and setting Zd = 1 



_Lp*(°) ( 



|q|) 



2^2 7 |q| 



4|q|7 - (3|q| 2 + 1QT 



tan 




At this order the deuteron quadrupole moment is (j,q 



(LO) 



Vsd 



(16) 



0.335 fm , deter- 



mined from the |q| — ► limit of eq. (|16|) (see also Ref. 

In the ^-parameterization at NLO for Fq corrections proportional to Zd — 1 appear. 
Also at NLO there is a contribution from a four-nucleon-one-quadrupole-photon local oper- 
ator whose coefficient is not related by gauge invariance to the A^ A" scattering phase shift 
data | 36fl . This coefficient must be determined from the deuteron quadrupole moment. An 



overall shift in the value of Fq, which we denote by 5/j.q, results from the inclusion of this 
operator. The total NLO contribution to Fq is 



J_,, 

m Q 



|q|) = SfJi Q +{Z d -l) 



^ (0) (|q|)-^ 



Ml 



{LO) 



(17) 



The numerical value of 5fJ,Q required to reproduce the measured value of the deuteron 
quadrupole moment is 5liq = —0.0492 fm 2 . The EFT analyses of Refs. suggest 
that a counterterm of this size is expected from renormalization group arguments and natu- 
ralness. As it is the contribution of a two-body quadrupole charge operator to Fq it simply 
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100 200 300 400 

k (MeV) 

FIG. 1. The deuteron quadrupole form factor. The dashed, dotted and solid curves correspond to 
the LO, NLO and ISpLO contributions in the z -parameterization. The dot-dashed curve corresponds 
to a calculation with the Bonn-B potential pzjj in the formulation of . 



does not appear in effective range theory. In Refs. [p^ f41j it was argued that the absence of 
such a term could be responsible for "modern" potential models' persistent underprediction 
of hq by ~ 5 % in impulse approximation |^ , |50| , |5ll . The preeminent role played by this 
counterterm, as compared to others in the NN effective field theory, is because /iq is more 
sensitive to short-distance physics than many other deuteron observables. 

As in the charge form factor, the only correction to this result at N 2 LO comes from the 
finite-size of the nucleon. 

^ (2) (|q|) = -~ rl |q| 2 ^a (0) (lql) , (is) 

All the nuclear effects have already been taken care of by the demand that Z d be replicated 
exactly. 

The form factor computed at LO, NLO and N 2 LO is shown in fig. [I], together with a 
potential-model calculation which used the Bonn-B potential in the formulation of . 
The difference between -Fg(|q| = 0) for the two higher-order EFT(-f) calculations and the 
potential model calculation shown in fig. [IJ, occurs because the calculation using the Bonn-B 
potential does not produce the correct deuteron quadrupole moment. 

We can also easily compute the deuteron quadrupole radius, defined by 



/ r 2\ = _ 6 ^k 



(19) 



At N 2 LO there is no counterterm contribution to (tq), and so we expect that it will be less 
sensitive to short-distance physics than \iq. Specifically: 

9 i] sd ( 80 2 2 



^ = ^( z < + f^<) - L40fa4 ■ (20 » 

At next order a contribution proportional to 7] 2 sd enters, and a two-body counterterm also 
appears. The error in our result for (r^) is dominated by the counterterm, which produces 
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a relative error (•y/m n ) 5 l/r] si i ~ 15 %. This is roughly the discrepancy between eq. ( p0|) and 
the results obtained from potential models with similar values of r] S( i. 

Similar results to those shown here for Fc and Fq can also be obtained for the deuteron 
magnetic form factor, Fm- It too has good convergence properties in the ^-parameterization 
of EFT (ft). 

We now leave the form factors for elastic electron-deuteron scattering and instead turn 
our attention to a quantity which is very much dominated by the long-distance physics of 
the deuteron. The electric polarizability of the deuteron has been computed in EFT (ft) up 
to N 3 LO in the p-parameterization [RRKRO: 



a 



(p) 

EO 



aM 



N 



32 7 4 



1 + Ipd 

0.377 + 0.153 
0.616 fm 3 , 



+ 



2 2 

7 Pd 



+ 



where Dp 



b) 



+ 2 



0.062 + 0.022 



SM 2 N 
+ 0.0006 



+ 



m n t 

3tt 
0.0036 



D, 



(21) 



■r & + ^§ ^C\ r2> = —1.51 fm a is the combination of P-wave 

coefficients that contributes to the El amplitude in np — > drj. The dominant corrections 
here always arise from the perturbative expansion of Z&. The relativistic correction, for- 
mally N 2 LO in the momentum expansion, gives a contribution which numerically is of much 
higher order. From now on we will ignore it. The uncertainty in this calculation in the 
p-parameterization is then set by the omission of 7 4 p^ and higher. This constitutes an 
uncertainty of ~ 2.5% due to the expansion of Z d , and leads to Aa^J ~ ±0.015 fm 3 . 
In contrast, the ^-parameterization gives 



a 



£0 



aM 



N 



32 7 4 



1 + (Z d 



+ + 



3tt 



0.3770 + 0.2605 + - 0.0036 



0.6339 far 



(22) 



rapid convergence indeed. The insertion of the P-wave operators also picks up factors of Z d 
when computed to higher orders and therefore we provocatively write 



a 



(*) 

B0 



aM, 



N 



32 7 4 

0.6314 fm 3 



1 + 



M N f 
3tt 



Dp + 



(23) 



where the ellipses denote higher-order terms, which include a N 4 LO contribution from S — D 
mixing proportional to rf sd . The first contribution from photon- four- nucleon vertices which 
are not constrained by NN phase shifts or gauge invariance also occurs at N 4 LO. We expect 
that these two effects combined may shift the value of oeo by ±0.0015 fm 3 . 

This result is commensurate with the good agreement between "modern" potential model 
calculations which give a^o — 0.6328 ±0.0017 fm 3 [54],|55|. The deuteron electric polarizabil- 
ity is very precisely predicted once the NN phase shifts are properly described. EFT thus 
provides a very accurate description of this quantity, as long as one demands that the tail 
of the deuteron wave function be well-reproduced at low order in the EFT expansion. This 
will be generally true for any low-energy elastic process on the deuteron. Examples of such 
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processes already computed in the p-parameterization include Compton scattering ||T6| , |T7 
and the reaction vd — > vd . 

We now turn our attention to reactions where the deuteron is only present in the final or 
initial state. In particular, we wish to see if the ^-parameterization improves the convergence 
of these reactions, in which only one of the asymptotic states involves a deuteron. 

The radiative capture process, np — » dj, has been studied in great detail in both EFT (ft) 
5E| , |3"8|J4"2] | and the theory with pions |SIJ . For incident cold or thermal neutrons, the rate for 



this process is dominated by Ml capture from the 1 Sq channel via the isovector magnetic mo- 
ment of the nucleon. The amplitudes for capture from other channels are strongly suppressed 
in this kinematic regime. Forthcoming experimental data from Grenoble has stimulated in- 
terest in the calculation of the isoscalar Ml and E2 matrix elements in np — > d'y using 
different formulations of effective field theory [[J5|,[8J. In Ref. [3SJ these amplitudes were 



computed in the p-parameterization of EFT^). We now wish to examine them in the z- 
parameterization. The isoscalar Ml amplitude has been computed up to NLO, and at that 
order it is the same in both the p and z parameterizations, due to the vanishing contribution 
of the one-body operator. However, the isoscalar E2 amplitude is affected by the choice of 



parameterization. Up to N LO the amplitude Xe2 s , as defined in [38], becomes 



AV2 10 V 8 V ; 32 v ; J 40 1 ' ' V ; 

when the ^-parameterization is used. Here E[ is found by fitting the shape of the S-D 
mixing parameter e x and turns out to be E^f 1 = —2.880 fm 4 [|38[|. 

Numerically, this leads to a series for the ratio of the isoscalar E2 matrix element to the 
experimental value of the dominant isovector Ml matrix element: 

-J^- = - [ 1.565 + 0.693 + 0.209 ] x 10" 4 = -2.47 x 10" 4 . (25) 
Xmi v 

While the central value is essentially unchanged from that found in Ref. [[55J we now estimate 
the uncertainty in this ratio to be ~ 6%, much smaller than the ~ 15% quoted there. The 
series in eq. ( P5|) appears to be converging faster than the analogous expression obtained 
using the p-parameterization ||38|| . The reasons for this are not entirely clear, since the naive 
expansion parameter, Z d — 1 = 0.690, is actually larger than 7p^ = 0.408. However, it seems 
that, in general, an expansion of \/l + — 1 in powers of Z& — 1 converges faster than one 
of (1 — 7Pd) -1 ^ 2 in powers of 7p^. It will be interesting to see how far the ^-parameterization 
improves the convergence of EFT calculations of the astrophysically-important reaction 
pp — > de + u e . Kong and Ravndal have pointed out that the main obstacle to agreement 
between EFT and potential model calculations of this reaction is the EFT's failure to exactly 
reproduce the overall factor of \fZ~2 which appears in the astrophysical ^-factor ||28|| . 

Finally, we turn our attention to iVW scattering. In this case we merely wish to reas- 
sure ourselves that the z-parameterization does no worse than the p-parameterization. We 
examine the 3 Si phase shift So up to N 2 LO. Taking the expression for the coefficients ( |10|) 
and computing the amplitude, and hence the phase shift, up to a given order in Q, we find, 
successively: 

5 = vr - tan" 1 f -) - (Z d -1)^- + (Z d -1) 2 ^- , (26) 
V7/ 2 7 4 7 
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1 00 200 300 

k (MeV) 

FIG. 2. The 3 S\ phase shift as a function of nucleon momentum |k|. The dashed, dotted and 
solid curves correspond to the LO, NLO and ISpLO contributions in the z -parameterization. The 
dot- dashed curve corresponds to the Nijmegen partial wave analysis 



where the first two pieces are order Q°, the second is order Q, and the third order Q 2 . 
Relativistic corrections are again omitted. This phase shift is plotted in fig. [| At N 2 LO it 
agrees very well with the Nijmegen partial wave analysis ||. The agreement is certainly as 
good as that obtained in the p-expansion ||36|| . It is very encouraging that even for momenta 
which are large in the context of this EFT, i.e. k ~ m^, the convergence of the Q-expansion 
for iVTV scattering is good if the ^-parameterization is used. Given this result, and that 
observables involving the deuteron appear to converge significantly faster, it is clear that 
the ^-parameterization of the coefficients in the Lagrange density (|5]) is better than the p- 
parameterization that has been used up until now. Its use may also improve the convergence 
of three- nucleon observables such as the quartet phase shifts p3|,33]. 

EFT(tF) is the most general theory in the two- nucleon sector consistent with the sym- 
metries of the strong and electroweak interactions in the very low-energy regime. As we 
explained above, the only change implemented here, vis a vis previous work on this theory, 
is to ensure that the tail of the deuteron wave function, in both the S and D-states, is 
correctly reproduced at low order in the EFT when elastic processes are calculated. That 
this improves the convergence of the EFT is not surprising. By choosing slightly different 
constraints to determine the coefficients in the Lagrange density we have been able to write 
observables in terms of the normalization constants and r] s d, rigorously defined in terms 
of the residue of the scattering amplitude at the deuteron pole. If the coefficients are fit 
in this way then the full impulse approximation result in effective range theory is obtained 
for elastic processes at NLO in the EFT — at least up to leading order in r) s d- NLO EFT(jc) 
calculations of elastic processes get the long-distance piece of the nuclear dynamics correct. 
However, EFT (ft) calculations then go beyond effective range theory, in that they systemat- 
ically include two-body currents which account for the short-distance physics not included 
in the asymptotic wave functions. The advantage of using the EFT is that these contribu- 
tions, which are not constrained by NN phase shift data, have sizes set by the underlying 
short-distance physics. The best example of this is the four-nucleon-one-quadrupole-photon 
operator which contributes at NLO to the quadrupole moment. 
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